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Abstract. - The one-dimensional Hubbard model at half-filling is studied in the framework 
of the Composite Operator Method using a static approximation. A solution characterized 
by strong antiferromagnetic correlations and a gap for any nonzero on-site interaction U is 
found. The corresponding ground-state energy, double occupancy and specific heat are in 
excellent agreement with those obtained within the Bethe Ansatz. These results show that the 
Composite Operator Method is an appropriate framework for the half-filled Hubbard chain and 
can be applied to evaluate properties, like the correlation functions, which cannot be obtained 
by means of the Bethe Ansatz, except for some limiting cases. 



Introduction. New materials whose physics is dominated by electron correlations in 
narrow energy bands are actually a challenge for solid state physicists. The treatment of such 
systems is not trivial due to the competition between itinerant and localized behaviour of the 
electrons in these bands. The Hubbard Hamiltonian Q is regarded as the simplest model 
which can give us the basic understanding of the effects of strong electronic correlations. In 
particular, its one-dimensional (ID) version is interesting for several reasons. On the one hand 
it is exactly integrable, by means of the Bethe Ansatz Q|. On the other, it could be applied to 
study real quasi-lD systems like the KCP and TTF — TCNQ salts |3j, the doped spin Peierls 
chains and the Cu-0 chains of the high-T c superconductors, whose underlying physics is 
directly related to the low-dimensionality of the system. 

By using the Bethe Ansatz and following a method developed by Yang Q, Lieb and Wu 
evaluated exactly some ground-state properties of the ID Hubbard model at half filling ||. 
Two years later, the spin and charge excitation spectra were obtained, also within the Bethe 
Ansatz, by Ovchinnikov M. For arbitrary electron density, the Bethe Ansatz coupled integral 
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equations for the charge and spin distribution functions cannot be analytically solved except 
in some limiting cases Q and numerical calculation is needed. In this way, Shiba || evaluated 
the ground state energy, local magnetic moment and magnetic susceptibility as a function 
of the electron density. The finite-temperature formalism of the Bethe Ansatz developed by 
Takahashi [jlO| provides the evaluation of the thermodynamic properties |Tl| . Further studies 
of the ID Hubbard model within the framework of the Bethe Ansatz have clarified for instance 
the thermodynamics in the presence of a magnetic field [[[2) . 

Therefore, we see that the Bethe Ansatz allows to exactly evaluate many quantities, thus 
providing a wide picture of the physics of the Hubbard chain. Nevertheless, this Ansatz cannot 
be regarded as a complete framework since many important properties, like the charge and 
spin correlation functions and the spectral properties, cannot be extracted from the exact 
Bethe Ansatz wave function except for some limiting cases (e.g. half filling, U — * oo, to = 0) 
[ jl3| . Otherwise we must address to other numerical H] or analytical jL{| approaches. 

The Bethe Ansatz is a very useful test for any approximation scheme, whose reliability 
can be checked by computing quantities that are given exactly by such Ansatz. In this 
sense, many of the available analytical methods do not give satisfactory results. On the 
other hand, all numerical techniques present some inherent problems, namely the small size 
of the clusters and the impossibility of reaching very low temperatures. Motivated by this, 
we have studied the ID Hubbard model at half filling by means of an analytic approach 
that has proved to be adequate for studying other strongly correlated models |l(| . The results 
obtained for the thermodynamic properties are in good agreement with the numerical data fin\ , 
also some anomalous thermodynamic and magnetic behaviours observed in high-T c cuprate 
superconductors have been successfully explained |i~7|, [ll). In this calculation scheme, called the 
Composite Operator Method (COM), the long-lived excitations of the system are described 
by an appropriate combination of the standard fermionic field operators. The properties 
of the new fermionic fields are self-consistently determined by the dynamics. To fix the 
internal parameters that appear, some symmetry requirements, like the Pauli principle and the 
particle-hole symmetry, are imposed. This procedure permits to recover symmetries that are 
badly violated by other approaches Jl9|], and thus is expected to provide a better description 
of strongly correlated systems |^0| . 

In this letter we analyze the half-filled infinite Hubbard chain within COM in the static 
approximation, where finite life-time effects are neglected. Although we are considering 
a paramagnetic (PM) ground state, we find a solution of the model which shows strong 
antiferromagnetic (AF) features. We present a detailed analysis of this AF-like solution. 
We calculate the energy and double occupancy of the ground-state, and the specific heat 
as a function of temperature, and compare them to the exact results obtained by the Bethe 
Ansatz. Excellent agreement is found. 

Method. - We consider the well-known Hubbard Hamiltonian: 

h = ^ ct (*) c Cj) + u "t (*) n i w - » E n W C 1 ) 

ij i i 

where (i) = ( cj (i) , cj is the electron operator on the site i in the spinor notation, n a (i) 

is the charge-density operator for the spin cr, /x is the chemical potential introduced to control 
the particle density n — {c i (i)c(i)) , and U is the on-site Coulomb interaction. Considering 
only nearest neighbours and taking the lattice constant as one, the hopping matrix for the 
chain is 

U,=-1t^Y. elk(l ' 3)co < k ) ( 2 ) 

k 
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In the case of the Hubbard model, a natural choice for the composite operators is the Hubbard 
doublet W = where 

£t(i) = c t(i)[l_„(i)] 

r,t(i) = (:*(») n(i) (3) 

This operators describe the hopping of an electron to an unoccupied and to an occupied site i, 
respectively. Considering a two-pole approximation | pl[ and a PM ground state, the Fourier 
transform of the single-particle retarded thermal Green's function may be written in COM as: 

^ lu - Ei (k) 

The energy bands are given by Ei(k) = R(k) + (— Q(k) where 

i?(fc) = u-2tcos(k) 




m(fc) = 2t[A + cos(fc) (p- h)] 

and/i, 72 are the diagonal matrix elements of the normalization matrix 7 = F.T. ^{\f r (z), ^(j)}). 
The spectral moments o"^(fc) that appear in the Green's function (H) are given in ref. pl[ . 

As shown above, the single-particle thermal Green's function (^) depends on the external 
parameters t, U , n and T (temperature), and three internal parameters /i, A and p. A and 
p are the following intersite charge correlation function and intersite charge, spin and pair 
correlation function, respectively jl6| po| . 

A = (t a (i)ZHi)) ~ (v a (i)vHi)) 

p = j(n^(i)nM-{(^(i)ci(i)rc\(i)c\{i)) (6) 

The superscript a indicates the field on the first neighbour sites and n M (i) = c) (i)tj ^c{i) is the 
charge- (/i = 0) and spin- (/i = 1,2,3) density operator, where is the Pauli representation 
of the 57/(2) symmetry group. These parameters roughly produce (see eq. (Js|)) a shift in the 
bands (A) and a bandwidth renormalization (p) . Very different results are obtained according 
to how these internal parameters are fixed p0|| . In COM they are determined by solving the 
following system of coupled self-consistent equations, 




(7) 

with S lS = (%ml(i)) and S« 5 = (*))• 

The first two equations are obtained from the existing relations with the elements of the 
Green's function, and the third one comes from requiring the satisfaction of the Pauli principle 
at the level of matrix elements (see ref. [|l6[ |o) for details). 

Once the internal parameters are determined, the evaluation of the physical quantities 
is straightforward. In this paper we study the single-particle band structure, the energy and 
double occupancy of the ground state, and the specific heat. The band structure is determined 
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Fig. 1. - Band structure of the half- filled Hubbard chain at T = for the indicated couplings. 



from eq. (j^)). The ground-state energy per site is calculated as the thermal average of the 
Hamiltonian and is given by 

E = 2t (g 1 - (1 - n)UF 1 + + ^U{1 - G - UF ) (8) 

where the functions F ni G n and B n (n = 0, 1, • • •) are defined in ref. pc{] . 
We use the thermodynamic relations E = F + TS, 

F(T, n) = £ M (T, n') dn' S(T, n) = - £ ^ dn' (9) 

with F the Helmholtz free energy and 5" the entropy, to determine the specific heat, which 
reads as 

As shown for the 2D case [0 , the temperature derivatives of the chemical potential can be 
expressed in terms of the internal parameters and calculated once the self-consistent equations 
(0) are solved. Finally, the double occupancy is calculated by deriving the free energy with 
respect to U. 

Results. - Considering a PM ground state, we find a fully self-consistent solution char- 
acterized by a negative value of the p parameter. The results for the energy spectrum are 
shown in fig. [I]. As we can see, this solution presents a typical AF band pattern; namely, 
a first excitation at k — ±-|, a very narrow bandwidth of order J = ^j- (the AF exchange 
interaction), and a quasi- halved Brillouin zone. Such an AF-like band structure is directly 
related to the negative sign of the p parameter, which is responsible of the general band shape 
[see eq. (§]. In then gure, the energy is measured with respect to the chemical potential. The 
solution exhibits a gap in the excitation spectrum for any non-zero value of U, in agreement 
with the Bethe Ansatz solution. The rate at which the gap increases coincides with Bethe 
Ansatz down to U «4. 
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Fig. 2. - (a) Ground-state energy and (b) double occupancy of the half-filled Hubbard chain in COM 
(solid line) and Bethe's Ansatz (dashed line). 




Fig. 3. - Specific heat versus temperature of the half-filled Hubbard chain in COM for U — 8. 



We calculate the ground-state energy and double occupancy as indicated in the previous 
section. These quantities, as a function of U, are shown in fig. || together with Bethe's Ansatz. 
The excellent agreement obtained is probably related to the opening of the gap mentioned 
above. Such a good agreement with Bethe Ansatz is not reached by other analytical approaches 
like the Gutzwiller approximations or the ladder and self-consistent ladder approximation 
[ fi"5| . In particular, these approaches fail to reproduce the correct asymptotic behaviour of the 
ground-state energy. As shown in panel (6), the double occupancy goes to zero as U — > oo. In a 
system of itinerant electrons with an average of one electron per site, the kinetic energy must 
also have this asymptotic behaviour, because any electron hopping would lead to a double 
occupied site. Therefore, electrons localize at infinite U, where the half-filled Hubbard chain 
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is equivalent to the spin-i AF Heisenberg chain that describes a system of localized spins. 

We calculate the specific heat of the system by means of eq. ([!(]) . The result as a function 
of temperature and for a ratio ^ = 2 (coupling/bandwidth) is shown in fig. ^. A double-peak 
structure appears, with a first peak located at T s» 0.25 and a second one at T « 1.6 (T is 
measured in units of t). This is in good agreement with the Bethe Ansatz calculation of Ref. 
[pl. The low-T feature is due to the spin excitations, since 0.25 is precisely the magnitude of 
the AF exchange parameter J = for U = 8 and a peak at such location is also found in the 
AF Heisenberg chain p2| . The high-T peak is associated to charge excitations. We can reach 
this interpretation by examining the chemical potential as a function of the electron density 
for several temperatures. At very low T, /i has a discontinuity at n = 1 of magnitude U that 
leads to the opening of a gap in the charge excitation spectrum, as we commented above. For 
increasing T the discontinuity in /j, is smoothed because the electrons can be excited across 
the gap. Therefore, charge excitations appear in the system at higher temperatures. Such a 
specific heat structure of the ID Hubbard model, with low- and high- T regions dominated by 
spin and charge excitations, respectively, is also obtained by numerical calculations on finite 



chains 23 . In the weak-intermediate-coupling regime < l) charge and spin excitations 



cannot be distinguished since they are of the same energy range, and hence the specific heat 
shows only one peak p3[| . 

Conclusions. Summarizing, we have studied the ID Hubbard model at half filling 
by means of the COM considering a PM state and a two-pole approximation. Within this 
approach we find a solution of the model that reproduces very well the exact solution given by 
the Bethe Ansatz. Namely, we obtain an insulating ground state characterized by strong AF 
correlations for any non-zero interaction U . The energy and double occupancy of this AF-like 
ground state is in excellent agreement with that of Bethe Ansatz. We have also calculated the 
temperature dependence of the specific heat of the system in the strong coupling regime. The 
locations of the spin- and charge-excitation peaks are consistent with the ones of the Bethe 
Ansatz. The approximation considered seems thus to be an adequate framework to study the 
physics of the half-filled Hubbard chain. It is of particular interest to apply this approach to 
get information about properties, like the correlation and spectral functions, which cannot be 
extracted from the Bethe Ansatz, except for some limits. 

The good agreement obtained for the ground-state energy and double occupancy within a 
static approximation is not surprising. As it is well known, such an approach gives a good 
description of the high-energy sector even though is unable to distinguish the low-energy 
features. Anyway, these features, are swallowed when integrating over the whole energy range. 
In this line of thinking, it is worth noticing the good description that we obtain for both the 
low- and high- energy sector of the specific heat. We suspect that this is due to our correct 
treatment of the Pauli principle at the level of matrix elements. In a physics dominated by 
strong electronic correlations, we believe that the satisfaction of the Pauli principle not only 
at the operator level, but also at the level of observation (relation among matrix elements) is 
crucial. 
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